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Abstract 

An approximate groundstate of the Anderson- Friedel impurity prob- 
lem is presented in a very compact form. It requires solely the opti- 
mization of two localized electron states and consists of four Slater 
states (Slater determinants). The resulting singlet ground state en- 
ergy lies far below the Anderson mean field solution and agrees well 
with the numerical results by Gunnarsson and Schoenhammer, who 
used an extensive 1/A^y-expansion for a spin 1/2 impurity with double 
occupancy of the impurity level. 

PACS: 85.20.Hr, 72.15.Rn 
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1 Introduction 



The properties of magnetic impurities in a metal is one of the most intensively 
studied problems in solid state physics. Although some of the experimental 
anomalies were already discovered in the 1930's it is still a subject of great 
interest. The work of Friedel PP and Anderson |2] laid the foundation to 
understand why some transition metal impurities form a magnetic local mo- 
ment while others don't. They considered a host with an s-band in which 
a transition metal atom is dissolved. The s-electrons can hop onto the d- 
impurity via the hopping matrix element Vgd- The ten-fold degeneracy of a 
real d-impurity is simplified and reduced to a two-fold degeneracy for spin up 
and spin down. If both states are occupied they repel each other due to the 
Coulomb exchange energy. This yields the Friedel- Anderson Hamiltonian 

N N 
HfA = { EiyCl^Cua + Eddldg + Vsdjl^) [dlCua + Cl^d„] ] + U Ud+Ud- (1) 

Here a finite s-band with N states is used. The c*^ and the d* are the 
creation operators of the (free) s-electrons and the d-impurity. The (io--states 
are assumed to be orthogonal to the s-states c*. 

Kondo |SI showed that multiple scattering of conduction electrons by a 
magnetic impurity yields a divergent contribution to the resistance in per- 
turbation theory. In the following three decades a large number of sophis- 
ticated methods were applied to better understand and solve the Kondo 
and Friedel-Anderson model, and it was shown that at zero temperature a 
Friedel-Anderson impurity is in a non-magnetic state. To name a few of 
these methods: scaling renormalization [H], Fermi-liquid theory [7], 
[H], slave-bosons (see for example P), large-spin limit [TU], [TTj, and the 
Bethe-ansatz ^21, QS]- For a review see ^3]. After decades of research exact 
solutions of the Kondo and Friedel-Anderson problems were derived |T2], J3] 
representing a magnificent theoretical achievement. 

The exact solution does not solve all questions. It uses a s-electron band 
with a constant density of states which extends from minus infinity to plus 
infinity (the cut-off is only performed at the end of the calculation). Further- 
more it is such a complex solution that only a limited number of parameters 
can be calculated and many non-critical or non-divergent contributions are 
neglected. For the majority of practical problems one uses approximate solu- 
tions. One particularly popular method is the large-spin method which will 
be discussed below. 
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While the single impurity problem is intensively studied and well under- 
stood the many-impurity problem and the periodic Anderson problem are 
still in a rather incomplete state ^Hl- Any simplified treatment of the single 
impurity may provide a new tool to improve the treatment of the latter. 

In this paper I wish to introduce a new method to treat the Friedel- 
Anderson problem. This approach can be best explained by a discussion 
of the ground state of the Friedel Hamiltonian. This is done in chapter 2. 
In chapter 3 the singlet state of the Friedel-Anderson model is derived. In 
chapter 4 the results are discussed. In the appendix some details of the 
calculation are summarized. 

2 The Friedel Resonance 

The author's method to treat the Friedel-Anderson problem can be best ex- 
plained by a discussion of the ground state of the Friedel (resonance) Hamil- 
tonian Hp^. For an s-band with N states and a d resonance Hp^ has the 
form: 

N N 

Hp, = e.clc, + EddTd + ^ Vsd{v) [d*c, + eld] (2) 

Since Hp, is identical for spin up and down, I will ignore the spin at the 
moment. 

As shown in ref. ^7], ^H] the exact ground state of Hp, with n (spinless) 
electrons can be written in the form 

n-l 

^Pr = [A'al + JJ a* $0 (3) 

i=l 

Here $o is the vacuum state and is a localized state which is built from 
the states of the s-band 

N 

«o = E (4) 

u=l 

The a*(l<i<A^ — 1) together with represent a new basis. The a* 
are orthogonal to Oq and to each other and their (A^ — 1) sub-matrix of the s- 
band Hamiltonian Hq = ^ eiyUi, is diagonal. (The construction of the states 
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{aQ,a*} is discussed in appendix A). The states a* are uniquely determined 
from the state Oq. Their form is 

N 

a* = J2(^l4 (5) 

In this new basis the (spin independent) Friedel Hamiltonian can be writ- 
ten as 

N-l N-l 

Hpr ='^E(i) a*ai + E (0) a^ao + ^ Vf^ (z) + a*ao] (6) 

i=l i=l 

N-l 

+ Edd*d + V^M[d*ao + aid] + ^ (0 + 

i=l 

where 

In the hamiltonian ® the first three terms represent the free electron 
hamiltonian. The Og-state represents an artificial Friedel resonance state. 
(APR state). It is interesting to note that the d*-siaXe and the localized 
OQ-state in ^ are on equal footing. The AFR state Oq is a sister state to the 
state d*. 

The terms with the matrix elements V^^ [i) and V^^ (i) yields the hopping 
between a* and d* and a* and Oq. For the state {A'a^ + B'd*) the individual 
hopping matrix elements cancel each other, making ^ the ground-state. 

Now I consider a conduction band and a d state with two spin components. 
Again I take the Friedel Hamiltonian as given by eq. (0). Then the ground 
state is given by the product of the spin up and down states of eq. This 
exact ground state can be written as 



^(o) = E.«^^.«^ 



fr 



E 



(7) 
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n-1 
cr,j=l 

n-1 

= [Aa^_aS+ + Bd*_a*Q^ + Cal_d\. + Dd*_d\\ J]^ a*^$o (8) 

iT,i=l 

where 

^+ + 5+ = 1 , A^^Bl ^ 1 

A = A+A_ , B = (9) 
C = A+B_ , D = B+B_ 

Each of the four states "if a, ^s, and ^'d is normahzed and they are all 
orthogonal to each other because they differ in the occupations of the d^ or 
dl-state. In Fig.l the four components ^a, and of the state 

are graphically shown. 
In the present case one has A^ = A_ and B^ = B_. With the right choice 
of Oq the above state is the exact ground state of the Friedel Hamiltonian. 
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Fig.l: The four Slater states used in the present calculation, 
ponent either the (i*-state or the AFR state is occupied. 



In each com- 



3 The Friedel- Anderson Model 

In the presence of a Coulomb exchange interaction U one can use the general 
form of the state "^ms as an approximate ground state. This state consists 
of four Slater states and possesses considerable flexibility. One may drop the 
condition = A_ and = B_ and use A, B, C, D as free parameters 
fulfilling the condition 

\Af + \Bf + \Cf + \Df = 1 

and optimizes these coefficients. Far more important is the optimization of 
the states ao+ and Oq-. With ao+, cio~ the full new bases , are 

uniquely determined. Details of the numerical optimization are discussed in 
the appendix. 
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In the following calculations a half filled s-subband is used, i.e. the num- 
ber of occupied states in each s-subband is n = The energy expectation 
value of the state \1/ms is calculated. For zero Coulomb exchange energy U 
the spin up and down bases are identical as are the coefficients B and C. 
With increasing U the spin up and down bases are shifted with respect to 
each other, the coefficients B and C differ and a magnetic moment develops. 
The size of this moment and a comparison with the mean field solution is 
discussed elsewhere jl9]. 

3.1 The singlet state 

The ground state of the Friedel- Anderson problem is a singlet state. From 
"^Ms one can construct a mirror state by exchanging spin up and down. 
Combining the two states yields then a singlet state which I denote as ^55. 
It is given by the following expression 

= (Ti) T ^MS (it) 

n— 1 n—1 

i=l i=l 

n—l n—1 

T [Aal_^a*_^^ + Bd*_^al^^ + Ca*_^rf;^ + I^rfl^rf^J JJ a*+^ JJ a*_^<l>o 

i=l i=l 

The sign =p is chosen so that one obtains a singlet state. This state is not 
normalized and the "i?" and "C" components are not orthogonal to each 
other. This introduces some additional terms in the ground-state energy. 
Furthermore the matrix elements between the states ms] and msi become 
determinants of single electron matrix elements. This is discussed in the 
appendix. 

For the numerical calculation an s-band with electron states c* is used 
(in the following I denote single electron states by their creation operator). 
A logarithmic energy scale is used, as introduced by Wilson j3] in his Kondo 
paper. One uses a finer and finer energy scale when approaching the Fermi 
energy sp = ^- A brief description of these electron states is given in the 
appendix. 

The ground-state energy of the singlet state is shown in Fig. 3 as a function 
of U with Ed = —U/2 and compared with the mean field ground state energy. 
Its energy clearly lies below the energy of the mean field state 
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Fig. 2: A comparison between the ground-state energies of the Anderson's 
mean field calculation and the singlet state "^ss- 



4 Discussion 

4.1 Comparison with the large A^j-expansion 

A number of approximate solutions have been suggested in the literature 
in which a localized electron state forms a singlet state with the magnetic 
impurity, see for example [201, |2I];|in|- They have been suggested for the 
Kondo problem and the Friedel- Anderson model. One particularly popular 
approximation is the large Nj expansion (see for example jTHI, ^l])- In 
the large Nf expansion one assumes that the impurity has a large total 
angular momentum Jf {Jf because this method is often used for f-impurities) . 
The "spin" has then a degeneracy of Nf = {2Jf + 1) of the total angular 
momentum states. In the limit of infinite Nf (the large Nf limit) this method 
yields an exact ground state. For smaller spin, in particular for spin 1/2, one 
generally performs an expansion in powers of 1/Nf. 



MFT 
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Gunnarsson and Schoenhammer (G&S) [221 applied the large Nf method 

to a finite Coulomb interaction and spin -, including double occupancy of the 

impurity level. They calculated the ground-state energy in different approx- 
imations. G&S give their energy parameters in units of [eV] . I denote their 
parameters with the index "GS". These parameters are: half the bandwidth 
Bgs, the d-state energy E^^gs, the Coulomb energy Ugs- For the s-d-hopping 
transition they use an elliptic form 



where g (e) is the density of states (per spin). All these parameters are ener- 
gies or potentials. By dividing these energy parameter by Bgs one obtains 
the appropriate parameters for the present calculation. When the numerical 
calculation is completed the resulting ground-state energy must be multiplied 
with Bgs for a comparison with GS's results. 

Gunnarsson and Schoenhammer calculated in ref. [22] the ground-state 
energy (for Nf = 2) with the following parameters: Bqs = 6 eV, Ubs = 5 
eV, Ed^BS = —2.5 eV. They performed two calculations, one for Vgs = leV 
and another for Vgs = 2eV. The results are shown in tables ^ and |21 The 
first column gives the electron states used in the calculation (for details see 
[221 the appendix). The second column gives the calculated ground- 
state energies. In the third, fourth and fifth columns the symbols /°, 
give the probabilities for a d-occupation of 0,1,2. The last column gives the 
power of the (l/A^/)-expansion. The last row gives the results of the present 
calculation for the singlet state. The values for /o,/i,/2 agree perfectly. 
Also the ground-state energies are quite close with Egs = — 0.245e\^ and 
£"0 = — 0.239e\^ of the present calculation. 



For Vgs = 2eV^ the ground-state energy of the present calculation lies 
even below the value of the 1/iVj-expansion, as shown in table |21 

The state "g" requires the variation of more than 15^ , i.e. more than 10^ 
amplitudes. In the present calculation the singlet state requires the variation 
of 2A^=60 amplitudes. Keeping this in mind, the resulting ground state of 
the present calculation is rather compact. 




GS 



(10) 
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states 


Eo [eV] 


/o 


h 


/2 


param. 


0)+a+b 


-0.108 


0.001 


0.974 


0.025 




+c+d+e 


-0.238 


0.031 


0.938 


0.031 




+f+g 


-0.245 


0.034 


0.931 


0.034 




singlet 
state 


-0.239 


0.035 


0.931 


0.034 





Table 1: A comparision between the numerical results by Gunnarsson and 

Schoenhammer and the author for the case of A^^^ = 2. The parameters, given 
in the units used by GS, are Bgs = 6 eV, E^fis = —2.5 eV, Vgs = 1- eV, 
Ugs = 5 eV. The first column gives the states included in the large spin- 
method, the second column gives the groundstate energy. The third, fourth 
and fifth columns give the weight of zero, single and double occupation of 
the d-states. The sixth column gives the number of optimized parameters 
(amplitudes) in this calculation. 



states 


Eo [eV] 


/o 


/i 


/2 


expans. 


0)-Fa+b 


-0.628 


0.141 


0.778 


0.081 




+c+d+e 


-1.126 


0.140 


0.745 


0.115 




+f+g 


-1.217 


0.137 


0.732 


0.132 




singlet 
state 


-1.234 


0.140 


0.722 


0.138 





Table 2: A comparision between the numerical results by Gunnarsson and 
Schoenhammer and the author for the case of Vgs = 2eV. Everything else 
is identical to table 1 
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4.2 Properties of the Artificial Resonance State 

The states Oq^ and aQ_ are of particular importance of the present treatment 
of the Priedel-Anderson impurity. They determine the rotation of the s- 
electron basis in Hilbert space and therefore the solution of the problem. 
We analyze the composition of Oq^. in terms of the original s-state energies 
e,^. As discussed above, a^.^ is composed of the original s-basis c* with the 
amplitudes q;[J_|_ 

* ■^-vAT * 

In Fig.3a,b the coefficients q;°^ and Q;°_of the states Og^and a*Q_ are plotted 
for the parameters: Ucou — — —0.5, IV^^I — 0.1 and the number of 

s-statcs is = 32. One recognizes that the amplitudes at large absolute 
energies is very different for spin up and down. They are almost mirror 
images. 




Fig.3a,b: The coefficients for the AFR states and aQ_. The Wilson 
spectrum is used (for the region on the left side of the arrow the numbers v 
corresponds to an energy of E^, = —l/T"). The energy of the s-electron c* 
is {Ej,-\- E^_i) /2. On the right side of the arrow one has the corresponding 
positive energies. 

1 1 2 

For the analysis at small energies we plot the occupation density |agj_ | / {E,^ + E^^i) 
as a function of u. In Fig.4a,b these densities are shown for N — 32 and 



11 



N — 48. In the latter case the energy interval next to the arrow (zero en- 
ergy) is 1/2^ times smaller than for the left plot. Obviously the sub-division 
at the Fermi energy is not yet sufficiently small at the left plot for A'^ = 32. 




Fig.4a,b: The low energy occupation in the state Oq^ for different sub-division 

of the energy close to the Fermi energy. In a) the smallest sub-division is 
A = 2-^6 ?^ 1. 5 X 10"^ and in b) it is A = 2"^^ ^6x 10"^ While in the left 
plot A is not yet small enough, one observes in the right plot the occupation 
of ao_,_ has saturated. 

While amplitudes and occupations for large energies were rather different 
for aQ_|_ and aQ_ the occupation at small energies is almost identical. This is 
shown in Fig. 5 where the occupations of a^^ and aQ_ are plotted in the same 
figure. At energies close to the Fermi energy the occupation of Oq^ and aQ_ 
are essentially identical. On a linear energy scale at small energies the plots 
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in Fig.4b and Fig.5 are essentially identical. 




Fig. 5: The occupation of the spin up and down AFR states for = 52. At 
small energies (close to the center arrow) the two are almost identical. (The 
difference between $0+ and $0- is of the order of the width of the curve.) 



The average occupation density of the states aQ_,_ and aQ_ is 1/2 since the 
band ranges form -1 to +1. Therefore a density of more than 100 is quite 
large. 

The AFR states have weight at small and large energies. The weight at 
large energies is responsible for the large "perturbative" part of the ground- 
state energy. The weight at small energies is responsible for the anomalous 
behavior at low temperatures, the Kondo effect. 



5 Conclusions 

In this paper I suggest a very compact approximate ground state for the 
Friedel- Anderson impurity. Its center piece are two artificial resonance states 



13 



^0+, Oq- ^'^^ ^Pi^ down s-electrons. These are combined with the 

d-clcctrons for spin up and down d*^,d*_ into two-electrons states of total 
spin zero, i.e. [AaQ_aQ^ + Bd*_aQ^ + CaQ_d*\_ + Dd*_d*^] . Then for each spin 
a new s-electron basis {a*_|_} is built. These two bases are completely de- 
termined by the AFR states. Finally the {n — 1) lowest states of the two 

n-l 

basis are occupied yielding the s-electron background Yl (^ia^o- The com- 

i = l,(T 

positions of the AFR states calculated by numerical variation 

which rotates the s-electron bases in Hilbert space. This ansatz is exact for 
a spin degeneracy Nf of " 1" and infinity. 

The properties of the singlet state are investigated. Its ground-state en- 
ergy and the occupations /o,/i,/2 of the d-states are in good agreement 
with the extensive calculations by Gunnarsson and Schoenhammer using the 
large iVj-expansion. However, while in the large A^j-expansion one has to go 
to a large basis of states to obtain a good ground-state energy the present 
solution is extremely compact. 

The spectral composition of the two AFR states is quite interesting. Their 
composition is quite different away from the Fermi energy. Close to the Fermi 
energy one finds a large peak in the occupation density which saturates only 
for very small energies. This low energy occupation is essentially identical 
for the spin up and down AFR state. 

The structure presented here of the ground state allows a number of 
variations. Instead of using just two bases {a^^, and {a^^, a*_} one can 
use four or eight bases, for each of the ipx = A, B, C, D) a different one 
for spin up and down. These solutions improve the ground-state energy only 
slightly for the singlet state. 

A detailed analysis of the present solution is planed. For example, the 
construction of the triplet state and the calculation of transport scattering 
by the impurity in this ground state are desirable. Above all, it is of interest 
whether an extension of the presented solution can contribute to the periodic 
Anderson impurity problem. 
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A Some Details about the Numerical Calcu- 
lations 



A.l Wilson's s-electron basis 

Wilson [S] in his Kondo paper considers an s-band with energy values ranging 
from —1 to 1. In the next step Wilson replaced the continuum of s-states 
by a discrete set of states. This is done on a logarithmic scale. The dis- 
crete energy values are 1, 1/A, 1/A^, etc and —1, — 1/A, — etc where 
A is a parameter larger than one. (In this paper A = 2 is chosen). These 
discrete points are used to define a sequence of intervals: the interval u 
(for z/<N/2) is = -1/2"-^ < e < -1/2" = (there are equivalent 
intervals for positive ^-values where u is replace by (A^ — u) but we discuss 
here only the negative energies). The new Wilson states c* are a superpo- 
sition of all states in the energy interval {C,u-i,C,u) and have an (averaged) 

energy (^^ + ^^_i)/2 = 

spectrum continues symmetrically for positive energies. 

The amplitude of the states at the origin are chosen to be (pu (0) = 
i^u — ^u-i) / = 1/2'^+^. A state which is homogeneously composed of 
all energies in the full band has than an amplitude of "1". 

Therefore this choice of s*-states yields a dependence of the s-d matrix 
element Vgd (i^) on the state u. The essential advantage of the Wilson basis 
is that it has an arbitrarily fine energy spacing at the Fermi energy. 

A. 2 Construction of the Basis QjQ ^ CL^ 

For the construction of the state and the rest of basis a* one starts with 
the s-band electrons {c*} which consists of states (for example Wilson's 
states). The (i*-state is ignored for the moment. 

• In step (1) one forms an normalized state Oq out of the s-states with: 

N 

< = Y.^l< (11) 

v=l 



I.e. 



3 
"4' 



3 

'8' 



1 



' 16' ■■' 2-2^^/2 ' 2-2^^/2 



. This 
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The coefficients can be at first arbitrary. One reasonable choice is 



1/ViV 



In step (2) [N — 1) new basis states a* {1 < i < N — 1) are formed 
which are normahzed and orthogonal to each other and to a^. 

In step (3) the s-band Hamiltonian Hq is constructed in this new basis. 
One puts the state Cg at the top so that its matrix elements are Hqi 
and Hit). 

In step (4) the {N — l)-sub Hamiltonian which does not contain the 
state Og is diagonalized. The resulting Hamilton matrix for the s-band 
then has the form 



m 

Vfr{2) 



Vfril) Vfri2) 

E{1) 
E{2) 



\ VfriN-1] 











Vfr{N 






1)\ 



(12) 



E{N-1) J 



The creation operators of the new basis are given by a new set of {a*} , 
(0 < i < — 1). Again the a* can be expressed in term of the s-states; 
a* = c^l^l. After the state Oq is constructed the other states 

a* are uniquely determined. The additional s-d hopping Hamiltonian 
can be expressed in the terms of new basis and one obtains the Friedel 
Hamiltonian as given in eq. (jHl). The state "^ss is formed and its energy 
expectation value is calculated. 

In the final step (5) the state Oq is rotated in the A^ dimensional Hilbert 
space until one reaches the absolute minimum of the energy expectation 
value. In the example of the Friedel resonance Hamiltonian this energy 
agrees numerically with an accuracy of 10~^^ with the exact ground- 
state energy of the Friedel Hamiltonian |17j . 



A. 3 The effective s-d matrix element for the multi- 
electron states 

The calculation of the energy expectation value requires the calculation of 
many-electron matrix elements in different bases. We sketch here an example. 
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We consider the more general case that we have two wave functions — 
<-<+K=la<^o and = &o-0?;nri^L^o. Each is built from two 
different bases: {ao+, 0.4+ } , {o-o-) o-i-} and {60+) i {^o-i ^^i-} (only within 
this section the operators b^ b* are used for the AFR states to distinguish the 
different basis systems). The energy expectation value contains for example a 
matrix element of the form Here the s-d Hamiltonian H^^ can 

be expressed in any basis but for this matrix clement the representation 
is the optimal one. For the above matrix element one needs only the hopping 
for spin up (+) : 



7V-1 



i=0 



(13) 



The only term in H^^ which yields a non- vanishing contribution to (VJ^^) 

This matrix element contains 

(a) the multi-scalar product of the two n electron states for spin down 

F^"" = (K- nr=i' bi^^oK. nr=i' «:-*o) and 

(b) the matrix element Mff = (d; nr=7 b^^l Eio' Kf (0 nr=i' <+^o) 
The multi-scalar product is a determinant of order n containing the single 
electron scalar products between all occupied states. 



(K-K- > {KM- > 



■^(n-l)-l"0 



bl 



(n-l)-l'^l 



h* \n* 



(14) 



When the two AFR states are identical then the underlying matrix be- 
comes the unity matrix. 
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Part (b) yields 



Kr(o) 



in - 1) 

^l+l'^(n-l)+ 



(15) 



/)* Id* 

"(n-l)+l"0+ 



^(n-l)+l'^l+ 



^(n-l)+ I «(«-!)+ ^ 



A. 4 Details of the comparison with Gunnarsson and 
Schoenhammer's numerical evaluation 

Gunnarsson and Schoenhammer (GS) [221 applied the large A^/ method to 
finite Coulomb interaction and spin 1/2. They calculated the ground-state 
energy in different approximations. Since it is interesting to compare their 
results with the present calculation I sketch briefly the different states they 
considered. The corresponding graphical sketch of these states can be found 
in flgure 1 of ref. [22|- These states are collected in table ??. The flrst column 
gives GS's code for the state, the second column shows in which power of the 
1/A^/ expansion the state occurs, the third column gives the occupation of 
the d-level in the considered state, the fourth column the number of holes and 
electrons (above the Fermi energy) in the s-band, and flnally the flfth column 
gives the number of amplitudes (parameters) which one has to optimize in the 
numerical evaluation (again is the number of band states in the numerical 
evaluation). As an example the state "d" is part of the (lA'^/)^ expansion, 
it has, for example, the dj-state occupied, the s-band has one hole in the 
s|-band, another electron-hole pair is either in the sj- or s|-band. The total 
multiplicity of the state "d" is therefore 2 * (A^/2) * 2 (A^/2)^ The prefactor 
2*2 is replaced by "a" in column 5 (a > 1). 
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name 
of 

state 


power 
of 

expans. 


occup. of 
d-states 


s-band 


no of 
par am. 


|0) 




empty d-states 


half occupied band for spin | and J, 





a 




1 df or d| 


one f- or |-hole in 0) 


N/2 


b 




2 d-states 


f-hole and |-hole in 0) 


a {N/2Y 


c 




empty d-states 


one |- or | electron-hole pair in 0) 


a {N/2f 


d 




1 dt or d| 


two holes in 0) and one electron 


a {N/2f 


6 




2 d-states 


three holes in |0) and one electron 


a {N/2y 


f 


imr 


empty d-states 


two t- or I electron-hole pairs in 0) 


a {N/2Y 


g 




1 d| or dj. 


three holes in 0) and two electrons 


a {N/2f 


h 




2 d-states 


four holes in 0) and two electrons 


a {N/2f 



Table 3: Gunnarsson and Schoenhammer states 



G&S use a (different) exponential energy mesh of the form = ± [a — exp (xj)] . 
They use the value a — 0.2, and xi lies in the range (In {a) , In (a -|- Bos)) 
= (—1. 609 4, 1. 824 5) .This means Xi takes the values Xi = In (a) +i/N * 
^\i{a + Bgs) — In (a)]. G&S used for the values 9, 19, 29 and extrap- 
olated to — > oo. For this comparison I use the corresponding energy 
mesh and extrapolation. The only difference is that in my calculation the £j 
yield the energy frame and the energy states lie in the center between two 
Si whereas G&S used the as their energy states. After the extrapolation 
towards N ^ oo this difference should be negligible. 

The energy dependent s-d-matrix element V (e) adds a complication in 
the numerical evaluation. It varies strongly with energy. I average [V (e)]^ 
over each energy range. 
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